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Abstract 

A derivation of the Fluctuation-Dissipation Theorem for the microcanonical ensemble is pre- 
sented using linear response theory. The theorem is stated as a relation between the frequency 
spectra of the symmetric correlation and response functions. When the system is not in the ther- 
modynamic limit, this result can be viewed as an extension of the fluctuation-dissipation relations 
to a situation where dynamical fluctuations determine the response. Therefore, the relation pre- 
sented here between equilibrium fluctuations and response can have a very different physical nature 
from the usual one in the canonical ensemble. These considerations imply that the Fluctuation- 
Dissipation Theorem is not restricted to the context of the canonical ensemble, where it is usually 
derived. Dispersion relations and sum rules are also obtained and discussed in the present case. Al- 
though analogous to the Kramers-Kronig relations, they are not related to the frequency spectrum 
but to the energy dependence of the response function. 
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I. INTRODUCTION 



The relation between the fluctuations occurring in a system at equilibrium and dissipation 
effects dates back to Einstein [lj and his theory on Brownian motion. After that, Nyquist 
[2] derived a relation between the electrical resistance and voltage fluctuations in linear 
electrical systems. It was realized then by Callen and Welton [3] that such a relation could 
be proven for general linear dissipative systems using quantum mechanics. At that moment, 
the intuition of the authors, as described in the last paragraph of their Introduction, was 
that the relationship between equilibrium fluctuations and irreversibility would provide a 
method for a general approach to a theory of irreversibility and, indeed, this was the way 
pursued by Kubo [4] to achieve the theory of linear response. It is well established now 
that linear response theory gives a general proof of the Fluctuation-Dissipation Theorem 
(FDT) which states that the linear response of a given system to an external perturbation 
is expressed in terms of the fluctuation properties of the system in thermal equilibrium. 

Because of this deep relation between the FDT and linear response theory, it is worth 
noting that the response, as formulated by that theory, is given for any equilibrium ensem- 
ble. In other words, the response function can, in principle, be known not only when the 
system is initially in thermal equilibrium but also in another equilibrium state such as, for 
example, the microcanonical one. Therefore, the theory is quite general in the sense that 
the linear response of a system and its equilibrium fluctuations could be related to each 
other for any kind of equilibrium conditions. Indeed, fluctuation-response relations have 
been derived even in the context of stochastic systems 5), (| and non-Hamiltonian deter- 
ministic systems pi using linear theory. Perhaps the very first work concerning different 
equilibrium conditions from the thermal one in Hamiltonian systems is Ref.|8j, where the 
author shows that Kubo's formula can also be derived in the classical microcanonical en- 
semble as long as the thermodynamic limit is considered. However, for many and different 
reasons, much more attention was given for the statistical mechanics in the canonical ensem- 
ble than in the microcanonical one and the generality of linear response theory concerning 
different equilibrium conditions was not much explored. Of course, one could argue that the 
equivalence of the ensembles in the thermodynamic limit would be the reason for focusing 
just in the canonical ensemble, but recent developments have shown that there are indeed 
strong motivations to consider different equilibrium situations. For example, a path integral 
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representation for the quantum microcanonical ensemble [9j presented a few years ago was 
motivated by situations where the microcanonical approach may be more appropriate as for 
the description of systems at low temperatures or with a finite number of particles. The 
microcanonical ensemble has also been considered in relations between fluctuation and dis- 
sipation in systems far from equilibrium like the Crooks relation, where its microcanonical 
version helps to understand the connection between various of those fluctuation theorems 
[lo| . In Ref . fll| . a derivation of a microcanonical quantum fluctuation theorem was pre- 
sented. Considering the work performed by a classical force on a quantum system when it is 
initially prepared in the microcanonical state, the authors provide a relation that could be 
accessible experimentally to measure entropies. In the context of nanosystems, where the 
number of degrees of freedom constituting the environment is not always large enough, the 
microcanonical ensemble has also been considered. In Ref. Ijj], a quantum master equation 
was derived describing the dynamics of a subsystem weakly coupled to an environment of 
finite heat capacity and initially described by a microcanonical distribution. Finally, an 
analysis in the microcanonical state has also contributed to the recent debate about the 

n 

foundations of the canonical formalism |13| . 

The microcanonical ensemble implies a description of an isolated system. Therefore, 
one might ask how a relation between fluctuations and dissipation can be possible in a 
situation where no energy can be dissipated. In the present work, our goal is to explore 
the relation between fluctuations and response in microcanonical equilibrium conditions 
through the framework of linear response theory. As will be explained later, mainly after 
the development linear response theory, the name Fluctuation-Dissipation Theorem was 
associated with some relations which are analogs of the results presented here in the context 
of the microcanonical ensemble. That is the reason we took the freedom to call them also a 
FDT even in a situation where there is no physical mechanism for dissipation. The paper 
is organized as follows. In Sec. II the derivation of a FDT using linear response theory is 
presented and its validity is verified in a simple example. In Sec. Ill different dispersion 
relations and sum rules are derived in analogy with the usual Kramers-Kronig ones and 
their meaning is discussed. They are different because they are not derived in the frequency 
space, like the usual ones. Conclusions are presented finally in Sec. IV. 
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II. DERIVATION OF THE FLUCTUATION-DISSIPATION THEOREM 



We start by considering a system whose dynamics is given by a Hamiltonian H. An 
external force K(t) is applied to this system such that H is now perturbed by an external 
potential given by —AK{t). Following J4], the reponse function of the system due to the 
external force measured through an observable B is given, in linear response, by 



<p BA (X,t-t') = Tr ( p e {X) — \ A(0), B(t — if) 



Tr (p e (X)-[A(t'),B(t) 



(1) 



where [, ] is the commutator and p e (X) is the equilibrium density operator as a function of 
a macroscopic parameter A. One can also define the following correlation function between 
A and B 



(2) 



C BA (X,t-t') = Tr(p e (X)-{A(0),B(t-t')} 



Tr (p e (X)-{A(t'),B(t)} 



where { , } is the anticommutator. This function gives the spectrum of equilibrium fluc- 
tuations when the system is unperturbed. For the canonical ensemble, p e (X) = p e (/3) = 
e -pB jz{p)^ where (3 = (ksT)~ l , and the FDT establishes a relation between the spectra of 
4>ba and Cba- That means a relation between an equilibrium and a nonequilibrium quantity. 

Our goal here is to show that there is also a relation between ipBA and Cba in the 
microcanonical ensemble. First of all, let us start with the expression for the microcanonical 
density operator p e (A = E). Following [9(, we take it as 

5{E - H) 



Pe(E) 



Z(E) 



(3) 



where Z(E) = TrS(E-H). 

To derive the FDT, it is necessary to introduce an appropriate representation of 5(E — H) 
like, for example 

1 /"7+ioo _ _ 

' • (4) 



S(E - H) = — dzexp (E - H)z 

Z7TZ J'y—ioo 



Expressions ([T]) and ([2]) can be written now in the following way 



•>BA(E,t-lf) 



Z{E) 



Tr 



2m 



7+ioo 



7—400 



dze^ E - & > 



ih 



(5) 



(6) 



Z(E) " V 2m . 

It is important to note that, since the integrals in the complex plane are always convergent, 
the trace and integral signs can be interchanged. Doing that, it is convenient to define 
the following new quantities: (Pba{E, t — t') = Z(E)<pBA{E,t — t') and CBA{E,t — t') = 
Z{E)C BA {E,t - If) to obtain 

/"y+ioD 



^772 J y—ioo 

Cba(E, t-t') = — dze Ez E BA (z, t-t'), 

Z7T1 J •y—ioo 



(7) 
(8) 



where 



= Tr 



i(*0,B(f) 



A ,{A(f),5(f)} 



(9) 
(10) 



Since c/?^ and Cba are given as inverse Laplace transforms of xba and Fba, they also satisfy 
the following relations 



Xba(z,t) = / dEe- Ez <p BA (E,T), 
Jo 

Fba(z,t)= / dEe- Ez C BA (E,r), 
Jo 

where r = t — t'. We introduce now the Fourier transform of xba and Fba, 

I POO 

Xba(z,uj) = — dre iujt xba(z,t), 

Z7T J-oo 

F BA (z,u) = 7^1°° dre- WT F B A(z,T), 
and also the auxiliary function 

S A b(z,t) =Tr(e- Az A(t')B(t)). 



(11) 
(12) 



(13) 
(14) 



(15) 



Noticing that e Hz A(t') = A(t' + %zK)e Hz and using the cyclic property of the trace, we 
obtain 



Tr 



e- Hz B(t)A(t' + izk) =Tr e' Hz A(t')B(t) 



(16) 
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Using 



l r°° 

— / dr e~ tuJT Tr 



e- Hz B{t)A{t' + izTi) 



— / dr' e- iu>T> Tr 
2n J~oo 



e- Hz B(t)A(t") e zhw ,(17) 



zhw 



where t" = t' + izh and t' = t — t" , one obtains from ( TT5l) and ( TT6 



zTio; 



where 

S BA (z,u) = — / rfr'e-^ Tr 
Using (TTHl) in the Fourier transforms of (|T3j) and (|T4j) yields 

1 



(18) 



(19) 



Xba{z,u) = — S A b(z,u) - S B a(z,uj 
in L 

F BA (z,u) = - S ab {z,uj) + S ba {z,uj) 
Finally, from ( 1201) and ( I2TI) . we obtain 

Fba(z,u/ 



= Sba{z,u) 
S BA (z,u) 



(e zhuJ - l) 

' ih 

(e zhu) + l) 



(20) 
(21) 



,h . ( zhu\ _ . . 

^2 coth I — J XBA ^ Z ' U ' 



which is our quantum FDT. In the classical limit h — > 0, we obtain 

Fba{ z i w ) = — XBiifow), 



(22) 



(23) 



which is our classical FDT. One easily realizes from (|22|) and ( |23|) that the replacement of 2 
by /3 in those equations leads precisely to the quantum and classical versions of the FDT in 
the canonical ensemble. However, the physical nature of (I22p and (1231) can be quite different 
from that in the canonical case. Let us consider, for example, in the classical regime an 
ergodic and small system, small in the sense that it is not in the thermodynamic limit. 
Then the micro canonical ensemble averages in ([1]) and ([2]) can be replaced by time averages 
whose behaviors are given by the dynamics of the system. Therefore, the fluctuations in this 
case happen due to the dynamics of the concerned system itself and not due to the coupling 
to a thermostat as in the canonical ensemble. From this point of view, it is surprising that 
there is a simple relation between the FDT in the canonical and microcanonical ensembles. 
Indeed, if one wants to compare both cases, the inverse Laplace transform in z should be 
performed on fl22|) and fl23|) since the canonical FDT consists of a relation between 4>ba(P, 



and Cba{P,w), keeping the original macroscopic parameter (3. For the classical case, this 
can be easily done using fl23l) . leading to 

Cba(E, uj) = - f E dE'(p BA (E\u). (24) 
to Jo 

For the quantum regime, the inverse Laplace transform should be performed on ff22l . It is 
not hard to imagine how different the result will also be from the canonical case. 

In addition to the pure meaning of the relation between response and fluctuations, one 
may wonder whether (1221) and (T23]) can be useful or not. We would say they can be useful 
in situations where the microcanonical ensemble can be applied and the thermodynamic 
limit is not satisfied. However, what we mean by usefulness is the possibility of applying 
the FDT in a context very different from the ones considered so far to obtain response 
functions from correlation functions and vice versa. If by useful one meant to go further 
and speak about, e.g., transport coefficients, then one would have to discuss more carefully 
the linear response theory in the microcanonical ensemble, especially because van Kampen's 



objections 14| can be more trickier in this case. The first objection, concerning the validity 
of the linearization, could still be answered as usual, we believe, by the argument of the 
stability of the distribution functions [3, [llj]. The second objection, concerning the origin 
of the decay of correlation functions which lead to finite transport coefficients, cannot be 
answered as is done sometimes in the context of the canonical ensemble by coupling to 



an environment 
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171 ] . The reason is simple: to use the microcanonical ensemble one 



assumes an isolated system. A possible answer in this case would be the instability of the 
dynamics jisl. \l9\. However, the question of what "dissipation" would mean in the present 
context of the microcanonical ensemble would remain. This is because, originally, the name 
Fluctuation-Dissipation Theorem comes from the fact that part of the Fourier transform of 
the response function is related to the power dissipated by the system when a time-periodic 
perturbation is applied to it. But for an isolated system there will be no dissipated power. 
On the other hand, the Fluctuation-Dissipation Theorem, mainly after linear response theory 
was developed, has been associated with an equation relating the frequency spectra of the 
response function and of the corresponding symmetric correlation function. In this sense, 
( 1221) and (1231) are analogous to Eq. (6.16) of Ref.J^] for the microcanonical ensemble and 



therefore we took the freedom of calling them Fluctuation-Dissipation Theorems as well. 
Although beyond the scope of the present work, a general and deep discussion of the subtle 



points mentioned above as well as of the linear response theory for the microcanonical 
ensemble would be of great interest and value. 



A. Example: the Harmonic Oscillator 

As an example, we would like to check (122]) and (1231) for a simple system whose response 
and correlation functions are known directly. In order to do that, we choose a simple 
harmonic oscillator. We consider the case A = B = X where X is the position operator. 
To perform first the calculation in the classical regime, we define the classical analogs of fl5]) 
and (jEJ) as 

tp(E, t-t')= J dx d Po 5{E- H(x , Pa)) {x(t'), x(t)}„, (25) 
C{E,t-t')= J dx dp 5{E - H(x ,p )) x{t)x{t'), (26) 

where { , } is the Poisson bracket with respect to the initial conditions (x ,p ) and x{t) is 
the solution of the classical equations of motion for the position. The averages above can 
be easily performed, leading to 

tp(E, t) = sin (u t), (27) 

muni; 

C{E,t) = ^ cos (lu t). (28) 



We can now calculate 



X(z,u) = — / dre-^ / rffie- B >(£,r), (29) 

Z7T J-oo JO 
1 roo roo 

F(z,cu) = — dre- iUT dE e~ Ez C(E,r). (30) 



2tt J-oo Jo 

The results are 



x(z,uj) = -iuz ~g(u) dEe~ Ez E (31) 



mu^ ' ' Jo 



F(z,uj) = -g(u) / dE e~ Ez E, (32) 



muj^ ' ' Jo 



where g(ui) = (l/2n) J^dre WT cos (u t). Therefore, 



/ 



F(z,uj) = — x{z (33) 



ZLU 



which agrees with (1231) . 



Quantum mechanically, we can calculate directly ([9]) and (fTUj) for the harmonic oscillator 
using the energy eigenbasis 

v-^ p ^sin (cut) . 

X (Z,T )=Y, e — > 34 

F( v )^ e *E,^! (35) 

where i£ n are the energy eigenvalues. Therefore, for the Fourier transform x(z, tu) we obtain 

w) = £ e~ £ " 2 ^ P("o + w) - %o - w)] • (36) 

Using (1221) and (I36I) . we obtain an expression for F(z,u). Inverting the Fourier transform, 
we get 

^/ \ f *h , f zhu \ cos (uj t) . . 

Since 

Ve^ = ^coth(^) -L - , (38) 
Equation (|37|) can be written as 

F{,,r) = -£e-^K^P- (39) 

which agrees with fl35|) . This verification of (T22]) for the quantum harmonic oscillator is the 
same as in the canonical ensemble case if z is replaced by (3. However, here (122!) still has to 
be transformed back to energy. 



III. DISPERSION RELATIONS AND SUM RULES 

In the canonical ensemble, it is possible to derive relations 
parts of the Fourier transform of the response function [l5, [2^]. Those are the so-called 
Kramers-Kronig relations and mainly they express a causality property contained in the 
response function. In the present case, dispersion relations also hold in the z-space because 
4>ba and Cba are defined for positive values of energy. Equations (TTTT) and (TT2"]1 imply that 
Xba(z,t) and F B a(z,t) are analytic functions in the half plane Re(z) > 7, where 7 is 
positive. Therefore, in this region 

( \ 1 I a Xba(z,t) 

XBA{z °> T) = 2^t dz ^^r- (40) 



aetween the real and imaginary 



Since hin^^oo \xba{z, t)\ = 0, we can close the integration contour with a semicircle in the 
half plane where Xba(z,t) is analytic and a line along Ke(z) = 7 and send the radius to 
infinity to obtain from (I40p the relation 



XbaKVoiT) = — P / dy 

m J-00 y -y 



(41) 



where the choices 2 = 7 + % and z D = 7 + iy Q were made. The right-hand side denotes 
the principal value of the integral. Writing xba in terms of its real and imaginary parts, 
Xba = x'ba + ^x"bAi Eq. (l4T]) leads to the following dispersion relations 



1 f°° 
Xba(Vo, t) = -P dy 

7T 



Xba(Vo,t) 



71 



00 
00 



dy 



x"ba{v, t ) 

y-y 
y-y a 



(42) 
(43) 



As it is usually done 201 ]. from the two relations above, it is possible to derive the moment 
sum rules, which, in this case, are related to the energy dependence instead of the frequency 
spectrum. The derivation of such sum rules is sketched in the Appendix. The results for the 
first three moments are shown below, where the subscript BA was dropped for convenience: 



<f(0,r) 
^)(0,r) 

^ (2) (0,r) 



7T J-00 
1 t°° 



7T J-00 

1 r°° 

TV J-00 



dyy 
dyy 2 



x"(y,r) + 
x'(y,r) + 



y 

y (1) (0,r) 

y 



where 



^ (n) (0,r) 



& 



dE n 



E=0 



(44) 
(45) 

(46) 
(47) 



The moment sum rules above are related to the asymptotic expansion of Xba with respect 
to z (which means low-energy behavior). For small values of z (i.e. high-energy behavior), 
one obtains the following sum rules: 



y? M) (o,r) 

^- 2 )(0,r) 
^ M) (0,r) 



1 

7T J-00 



dy 



X"(0, r) 

y 



dy- 2 [x'M + ^-^r) 

7r J-00 y z 1 

dy\ fx"(0,r) + ^ 2) (0,r) 



7T J-00 



(48) 

(49) 

(50) 
(51) 
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where 



roo roo roo 

(0,r)=/ dE x dE 2 ---\ dE n ip(E,r). (52) 

J J JE n -i 



The procedure shown in the Appendix can be repeated as long as the derivatives and 
the integrals (p(~ n ^ exist to derive higher-order moment sum rules. 

As for the sum rules in the frequency space, those above can be used to correct phe- 
nomenological expressions for yp(E,r). For example, if one assumes a functional form for 
the reponse function with some free parameters with respect to the energy dependence, one 
could determine them by imposing the sum rules for high- or low-energy behavior. The way 
to do that in the frequency space is shown, for example, in [3, [3]. Since the relations above 
are valid for any value of r, one could also have dropped the r dependence by setting r = 0. 
Then, it is easier to understand the meaning and the importance of the sum rules: the z 
spectrum of % ls given in terms of static quantities like ip^ n \0,r = 0) and <£>(~ n )(0,r = 0), 
which could be calculated quantum mechanically in terms of the commutation relations 



between A and B (see, for example, {4]). 

IV. CONCLUSIONS 

Using linear response theory, we presented a derivation of the Fluctuation-Dissipation 
Theorem in the microcanonical ensemble in both quantum and classical regimes. The the- 
orem is stated as a relation between the Laplace-Fourier transforms of the response and 
symmetric correlation functions. Although this relation is very similar to the one derived 
in the canonical ensemble context, it is valid, for example, in a situation where the fluctu- 
ations are very different from thermal ones, namely, fluctuations of an isolated system that 
is not in the thermodynamic limit. Therefore, the Fluctuation-Dissipation Theorem can be 
considered as a much more general relation and not constrained just to the context of the 
canonical ensemble. We believe this result can be very useful to calculate correlation func- 
tions from response functions (and vice versa) for systems in the microcanonical ensemble 

nn 

when they are not in the thermodynamic limit. In this sense, as mentioned in |a 112] (see 



also the references in 12j), the present work can be considered as an additional effort to 
apply statistical physics to small systems. Moment sum rules were also presented for the 
energy dependence and they could be useful to correct phenomenological expressions for the 
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response functions. 



APPENDIX: DERIVATION OF THE SUM RULES 

In this appendix we give a brief sketch of how to derive the sum rules presented in sec. 



III. For a careful derivation and deeper discussion about the subject, we refer to 20j. Our 
starting point is the function f(z) defined by 



POO 

f(z)= / dEe- Ez <P(E), (A.1) 
J o 



where z = 7 + iy is complex with its real part positive and 4>(E) is real. Therefore, f(z) is 

analytic in the half plane Ke(z) > 7 and it satisfies the following dispersion relations: 

1 r°° f"(v) 
f(y ) = -P dy±^L, (A.2) 
7T J -00 y -y Q 

f"(yo) = --P r dy^L, (A.3) 
ii J -00 y -y 

where f'(y ) and f"(y ) are the real and imaginary parts of f(y ), respectively. From (1A.2I) . 
we can write 

1 r°° f"(v) 

f(0) = -P dyLW (A.4) 
7T j-oo y 

and from flA.3j) multiplied by y D we obtain 

hm y f"(y ) = - T dyf'(y). (A.5) 
To calculate the left-hand side of (IA.5I) . we go back to flA.lj) and integrate by parts to obtain 



f(z) = ^ + - rdEe- E *^\E), (A.6) 
z z Jo 



where 



i<"»(0) = (^JiE) 



E=0 



(A.7) 



Therefore, 



and, from HAM and fjA.8j) 



hm zf(z)=<p(0), (A.8) 

\z\— >oo 



1 r 00 

0(0) = -/ dyf(y), (A.9) 

7T J-00 
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which is the first sum rule. To derive the next one, we define a new function 

g(z) = / dEe- E *(fP(E), (A.10) 
Jo 

which is analytic again for Ke(z) > 7. Therefore, g(z) obeys the same dispersion relations 
as f(z). Integrating (IA.10I) by parts yields 

0«(O) 1 



9W 



+ - 1" 'dEe- E *<t>®(E), 
z Jo 



from which we obtain lim\ z \-+ 0O zg(z) = 4> ^(0). By the same procedure as before, 

0«(O) = - / dyg'{y). 

7T J -00 



Since, from (IA.6[) and flA.lOp . 

roc 

g(z)= dEe- Ez ^ l \E) = zf{z)-m 
Jo 

and 

g'(y) = -yf"(y)-m, 

we obtain, from (IA.12I) and (IA.14I) . the second sum rule 



(A.ll) 



(A.12) 



(A. 13) 



(A. 14) 



(1) (O) 



dyy 



dyy 2 



f"(y) + 



0(0) 



y 



71 J-00 

Repeating the same procedure again, we obtain 

^(o) = -i 

71 J-OO 

and so on, as long as the (j)^(0) exist. 

A similar procedure can be applied to generate a different kind of sum rule [20j. Starting 
again with (lA.lj) . we integrate by parts in a different way now: 



f *( ^ (1) (o) 



y 



(A.15) 



(A.16) 



where 



From flA~2|) . 



f( z ) = -0(-i)(o)+ 2 / dEe- EK <fi- x \EY 

Jo 



roo roo roo 

0(-™)(O) = / dE x / dE 2 --- dE n <p(E). 
Jo JEi Je u -i 



n J-00 y 



(A.17) 



(A.18) 



(A.19) 
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and from flATTBl . f (0) = -^-^(O), so 

^(0) = 



7r J-oo y 



(A.20) 



We again repeat the procedure, as before, defining from (1A.20|) a new function g{z 



9{z) 



roo roc 

/ dE e~ Ez IE) = -<P { - 2) (0) + z dE e 
Jo Jo 



(E) 



m | ^(o) 



z z 

Since g(z) satisfies the same dispersion relations as f(z), we obtain 

Ti J-oo y 

Inserting the imaginary part of the second line of (1A.21D in (1A.22D leads to 



1 r°° 1 
( - 2) (0) = - / dy- f'(y) + ^- 1 \0) 



7r J-oo y 

Repeating the same procedure again, we obtain 



1 r°° 1 r 
0(-3)( O ) = - / dy-\f"(y) + y<j ) ^(0) 
n J-oo y* L 



and so on, as long as the jjk n -*(0) exist. 
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